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Electronic and magnetic properties of ribbon-shaped nanographite systems with zigzag and 
armchair edges in a magnetic field are investigated by using a tight binding model. One of the 
most remarkable features of these systems is the appearance of edge states, strongly localized near 
zigzag edges. The edge state in magnetic field, generating a rational fraction of the magnetic flux 
(4> — p/q) in each hexagonal plaquette of the graphite plane, behaves like a zero-field edge state 
with q internal degrees of freedom. The orbital diamagnetic susceptibility strongly depends on the 
edge shapes. The reason is found in the analysis of the ring currents, which are very sensitive to 
the lattice topology near the edge. Moreover, the orbital diamagnetic susceptibility is scaled as a 
function of the temperature, Fermi energy and ribbon width. Because the edge states lead to a 
sharp peak in the density of states at the Fermi level, the graphite ribbons with zigzag edges show 
Curie-like temperature dependence of the Pauli paramagnetic susceptibility. Hence, it is shown that 
the crossover from high-temperature diamagnetic to low-temperature paramagnetic behavior of the 
magnetic susceptibility of nanographite ribbons with zigzag edges. 
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I. INTRODUCTION 

Nanographites are nanometer-sized graphite fragments 
which represent a new class of a mesoscopic systems 
intermediate between aromatic molecules and extended 
graphite sheets. In these systems the boundary regions 
play an important role so that edge effects may influ- 
ence strongly the 7r-electron states near the Fermi en- 
ergy. A useful and simple system to investigate the 
electronic states of nanographites is provided by ribbon- 
shaped graphite sheets. The study of the electronic states 
of graphite ribbons based on the tight binding model re- 
veals that the edge shape - we distinguish between zigzag 
and armchair edges (see Fig.l) - leads to a strikingly dif- 
ferent properties of the states near the Fermi level. The 
ribbons with zigzag edges possess partly flat bands at the 
Fermi level corresponding to electronic states localized 
in the near vicinity of the edge. These localized states 
( "edge states" ) correspond to the non-bonding molecular 
orbital (NBMO) as can be seen examining the anal; 
solution for semi-infinite graphite with a zigzag edgel 
In contrast localized edge states and the corresponding 



flat bands are completely absent for ribbons with arm- 
chair edges. 

The localized edge states are of special interest in 
nanographite physics, because of their relatively large 
contribution to the density of states (DOS) at the Fermi 
energy. There is a tendency towards a Fermi surface 
instability which is important to examine. Previously, 
it was reported that based on the Su-Schrieffer-Hegger 
(SSH) model the electron-phonon interaction does not 
induce a lattice distortion because of the non-bonding 
character of the edge statesa. On the other hand, the 
electron-electron interaction on the level of an unre- 
stricted Hartree-Fock approximation (HFA) yields a fer- 
rimagnetic spin polarization pal the zigzag edges and an 
energy gap at the Fermi levelM. The effect of long-range 
Coulomb interaction on the edge state was examined us- 
ing the Par-Parier-Pople (PPP) model with the restricted 
HFA which does not allow any spin polarization. The 
conclusion was that long-range Coulomb interaction does 
not destroy the edge state as well as the flat bandsu. This 
results was further confirmed by first principle calcula- 
tions based on local density approximation (LDA)Q. 
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The stability of the edge states has been extensively 
investigated from various points of view. One remaining 
problem is the influence of an external magnetic field. 
How are the edge states affected? Is their NBMO charac- 
ter preserved? The answers to these questions will be im- 
portant for future studies of the electronic, magnetic and 
transport properties of nanographites. For this purpose 
we investigate here the magnetic properties, especially 
magnetic susceptibility, for the case of non-interacting 
electrons. 

The observed magnetic susceptibility \ is the sum of 
four components: (1) localized spin susceptibility Xspin, 
(2) diamagnetic susceptibility due to the core electrons 
Xcore, (3) Pauli paramagnetic susceptibility xp an d (4) 
orbital diamagnetic susceptibility Xorb. due to the cy- 
clotron motion of the itinerant electrons. 

Since we neglect electron-electron interaction through- 
out this paper, Xspin can be neglected. Furthermore, 
Xcore is unimportant for us, because it is small and basi- 
cally temperature independent. On the other hand, the 
Pauli paramagnetic susceptibility is related to the DOS 
at the Fermi level, which represents an important compo- 
nent in zigzag nanographite ribbons where an enhanced 
density of states appears at the Fermi level. Note that 
is negligible in armchair ribbons, aromatic molecules 
and graphite sheets, because their DOS is suppressed at 
the Fermi level. We will see below that the fact that 
the DOS introduced by the edge states is sharply peaked 
at the Fermi energy, xp introduces a very pronounced 
temperature dependence which is nearly Curie-like. The 
diamagnetic contribution to the susceptibility is very fa- 
miliar from the magnetic properties of graphite sheets. It 
is due to the orbital cyclotron motion of the electrons in a 
field with a finite component perpendicular to the plane. 
Naturally, this diamagnetic response is very anisotropic 
and only weakly temperature dependent. From this we 
can conclude that in nanographite ribbons with zigzag 
edges the susceptibility should consist mainly of these 
two competing contributions, xp an d Xorb.- Hence, a 
crossover occurs from a high-temperature diamagnetic 
to a low-temperature paramagnetic regime, where the 
characteristic temperature depends on the width of the 
ribbon and of the orientation of the external field. It is 
worth noting that the field direction is an important tool 
to distinguish the magnitude of the two components. 



II. ELECTRONIC STRUCTURE OF GRAPHITE 
RIBBONS IN MAGNETIC FIELD 



A. Harper Equation 



In this paper, we use a single-orbital nearest-neighbor 
tight binding model in order to study the electronic 
states of nanographite ribbons. This model has been 
successfully used for the calculation of electronic states of 



fullerene molecules, carbon nanotubes and other carbon- 
related materialso'S. The Hamiltonian is written as, 



H = ^t 

(■>■]) 



(2.1) 



where the operator c\ creates an electron on the site i, 
denotes the summation over the nearest neighbor 
sites. Here, we neglect spin indices for simplicity. In 
this model, the magnetic field B perpendicular to the 
graphite plane is incorporatecLin the transfer integral 
by means of the Peierls phased defined as 



(2.2) 



where <pij is given by the line integral of the vector po- 
tential from i-site to j-site, 



y = -/ 61- A. 



(2.3) 



The magnetic flux through the area S in units of the flux 
quantum (f> = — is 



aroundS 



.>(] 



ch 



(2.4) 



The structure of graphite ribbons with zigzag and arm- 
chair edges are shown in Fig. |l|, where we assume that all 
edge sites are terminated by hydrogen atoms. The rib- 
bon width iV is defined by the number of zigzag lines for 
the zigzag ribbon and by the number of dimer lines for 
the armchair ribbons. Since a hexagonal lattice can be 
divided by two sublattices, we call the A(B)-sublattice on 
the n-t/i zigzag or dimer line as nA (nB) site. We assume 
Landau gauge with A = (0, Bx, 0), where we define the 
translational invariant direction of each ribbon as the y- 
axis, and the x-axis lies perpendicular to y-axis. In this 
gauge, the unit cell of each ribbon could be taken as the 
rectangle shown in Fig. |l|. 




-V » VK • \vT ,VB 

FIG. 1. The structure of graphite ribbons with (a) zigzag 
edges and (b) armchair edges. The rectangle with the dashed 
line is the unit cell. 
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It should be noted that the same number N for both 
zigzag and armchair ribbons does not give the same rib- 
bon width, when the ribbons are measured by the same 
unit of length. Therefore, when we compare physical 
quantities of zigzag and armchair ribbons with a same 
width W, we will use the following definition 



W 



(JV-l)^o 



W z zigzag ribbons 
W a armchair ribbons 



(2.5) 



where a is the C-C bond len gth. 

Next, let us apply Eq.(2.1) to the graphite lattice and 
deri ve t he so called Harper equations. In order to apply 
Eq. ( |2.l| ) to the graphite lattice and simplify the formula- 
tion, we introduce the lattice transformation as shown in 
Fig. ||, which does not change the lattice topology. For 
convenience we will use this brick-type lattice structure 
when ever we perform real calculations. 




FIG. 2. The transformation from hexagonal lattice to brick 
type lattice, which does not change the lattice topology. 

The Peierls phase for the graphite ribbons is easily cal- 
culated from Eq.(2.3). The Peierls phase of graphite rib- 
bons with (a) zigzag and (b) armchair ribbons is shown 
in Fig. ||. For both cases, the Peierls phase is given by 
4>mB.nA = ^Tn(j)6 mn , where cj> is the magnetic flux through 
a plaquette in units of a quantum flux. 
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FIG. 3. The Peierls phase for Landau gauge on (a) zigzag 
ribbons and (b) armchair ribbons. 

Now, let us derive the Harper equation, which possess 
the translational symmetry along the zigzag axis. In or- 
der to apply Eq.(2.1) to zigzag ribbons, we define a new 
operator, c a (i), which create an electron on the site i in 
the unit cell a. The momentum representation of this 
operator in the y-direction is defined by 



c a {i) 



ikr c 



(2.6) 



where r a represents the position of the unit cell a. Here 
we also define a one-particle state 



m 

+ *mB(khl B (k))\0). 



(2.7) 



Inserting this one-particle state into Schrodinger equa- 
tion H\^) — e | , we can easily obtain the following four 
eigenvalue equations for the sites mB, mA and (m+l)A, 

e*mB = y {m+ i)A + fP****mA + e~ i2 ^^ mA , 

e^ mA = * (m _ 1)B + e i2 *^ mB + c- i2 ^^ mB , (2.; 

e*(m+iM = t roB + e i2ff f ** (m+1)B + e- i2 ^t>V im+1)B . 

Eliminating the A-sublattice sites, we obtain the differ- 
ence equation, 



^ m -x(k y ), (2.9) 

where A = e 2 — 3, a m (k y ) — 2cos(-y- + mncf)), b m (k y ) — 
2 cos(k y + 2mn<j>), and ^ mB was replaced by There- 
fore, our problem was reduced to a one-dimensional tight 
binding model with a superlattice potential of period 2q. 
Note that this equation does not include any boundary 
conditions yet. It may be applied to both graphite rib- 
bons and sheets by imposing the appropriate boundary 
conditions. In the following calculations, the factor rmrcj) 
will be replaced by (-^-^ — m + l) 7T<p to keep the energy 
band symmetric about k = for arbitrary magnetic flux. 
This replacements mean that the origin of the x-axis is 
set to the center of the ribbons. 

The spectrum is confined to values of A between —6 
and +6, i.e. —3 < e < +3. A close inspection shows 
that the following translations do not change the energy 
spectrum; e — > — e and <f> — ► (f> + n , where n is an arbitrary 
integer. For rational flux <fi = 2 , function with 

period 2q and b m is a function with period q. In addition, 
we must pay attention to the following symmetry of the 
energy bands in the Brillouin zone, 



' I k y H n 



2n 

9. 



e{k v ) 



(2.10) 



Here we should again note that these arguments do not 
depend on the boundary conditions. 

Similarly, we can also derive the Harper equations 
which includes the translational sy mmet ry along the 
armchair axis. In the same way as Eq.(2.9) the following 
equation is obtained, 



X^ m {ky) = * m+2 (fcj,) + a m ^ m+1 (k y ) 
+a m ^ m -i(k v ) + ^ m -2(ky) 



(2.11) 



3, a m (k y ) = 2 cos(^ + (to - §)7r<£)e -i ^. 
It is easy to c onfi rm that the same symmetry proper- 
ties as for Eq.(2J;) apply here too. Similarly, the factor 
a m will be replaced by a m {k y ) = 2 cos(-# + {^-^- — 2m + 
3)7T(/))e _17r ^ in order to obtain the symmetric energy band 
structure about k — for arbitrary magnetic flux. 
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B. Graphite Sheet 

In this subsection, we consider the electronic struc- 
tures of a graphite sheet in a magnetic field. This will 
afterwards become the basis to discuss the electronic 
structures of nanographite ribbons. As we have seen in 
the previous subsection, the tight binding model of the 
graphite lattice could be reduced to a tight binding model 
with superlattice potential of period 2q. In order to cal- 
culate the energy spectrum of the graphite sheet, we must 
treat the eigenvalue problem of a 2q x 2q matrix with the 
periodic boundary condition ^2q+i = e x2q ^>x, when the 
Brillouin zone is reduced to the magnetic Brillouin zone 



f q < k x < % and 



-7T < k y < 71". 



At the beginning, we consider the zero-field energy 
band of th e graphite sheet. By setting cf> — in Eqs.(^jj) 
and ( 2.11 ), the zero-field spectrum is easily obtained. 
From Eq.(p|), we find for the graphite sheet with trans- 
lational symmetry along zigzag axis, 



6£ = ±W3 + 2cos(k y ) + 4 cos 



cosOz). (2.12) 



Similarly, from Eq.( [2.11 ), the zero-field spectrum of 
graphites with translational symmetry along the arm- 
chair axis is 



= ±1 



2 008(2^) 



4cos ( -j- ) cos(fe x ). (2.13) 



The energy band structures for both Eqs.(2.12) and 
(2.13) are shown in Fig. ||, where k y is replaced by k and 
we have superposed all k x -values in the spectrum. We 
can find the degeneracy at e = in both figures, which 
originates from the K-poin^degeneracy of the band struc- 
ture of the graphite sheet! 





It -71 o It 

k k 
FIG. 4. The energy band structure of graphite sheet pro- 
jected to (a) zigzag axis and (b) armchair axis in the absence 
of magnetic field. 

The spectrum of the graphite sheet in a magnetic 
field is shown in Fig. [|, which was first calculated by 
RammalEj. We can easily find that the spectrum has the 
fine recursive structure of the Hofstadter butterfly. In the 
weak-magnetic field limit, we can clearly see the Landau 



levels. When the magnetic flux is getting larger, these 
levels form the Landau subbands because of the Harper 
broadening. As we pointed out in the previous subsec- 
tion, for the rational flux (f> — | , we can see 2q subbands 

with a reflection symmetry about e = and about <j> = \ . 
Interestingly, the degeneracy at e = exits for arbitrary 
flux, which confirms that the K-point degeneracy will not 
be destroyed. 




FIG. 5. The energy spectrum of the graphite sheet in a 
magnetic field. 
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FIG. 6. The energy band structure of graphite sheet pro- 
jected to (a) zigzag axis and (b) armchair axis for cf> — 1/4. 



As an example of the energy band structure in a strong 
magnetic field, we show the case of the graphite sheet 
projected to the (a) zigzag axis and (b) armchair axis for 
0=1/4 (Fig. [6]). In a strong magnetic field, the Landau 
levels change to 2 x 4 subbands due to Harper broad- 
ening and the effect of the lattice structure gets more 
important, so that each subband has the basic structure 
of the zero field case. Furthermore we find the symmetry 

2tt 
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C. Graphite Ribbon 



The energy band structures of graphite ribbons are ob- 
tained in terms of the Harper equation by imposing open 
boundary conditions. In the case of zigzag ribbons with 
width N, the boundary condition is ^n+i = = 0. 
However, we need a more careful t reat ment of the Harper 
equation at the edge site. In Eqs^^), it was not consid- 
ered that there would be no 0A(B) and (N+1)A(B) site. 
Including this fact, only for m = 1 and N, the Eq.( |2.9| ) 
has to be rewritten as 



+ {b m - l)V m (ky) 



+a m - X ^ m -x{k v ). (2.14) 




Theref ore we must replace b\ (6at) by b\ — 1 (b^ — 1) 
in Eq.(2.8) in order to include the condition that the OA 
and (N+1)A sites do not exist. Similarly, for the arm- 
chair ribbons, the boundary condition is \&./y+i = = 0, 
where we must replace b\ (6/v) by b\ — 1 (b^ — 1) in order 
to include the condition that the OA and (N+1)A sites 
do not exist. 



-71 71 -71 71 

k k 
FIG. 7. The energy band structures of zigzag ribbon with 

N = 50 for (a) <f> = 0, (b) cj} = £5 , (c) <j> = -L and (d) cf> = \ 

for (j> = i. The dashed line in (b) is the energy structure at 

= for comparison. 



We show the energy band structures of the zigzag rib- 
bon with N = 50 for = 0, ^, ^, \ in Fig. % (a)-(d). 
For = 0, the profile of band structure has almost the 
same structure as in the case of the graphite sheet as 
shown in Fig. [| (a). However, we can see partly flat 
bands at E = 0, which do not appear in the energy band 
of the graphite sheet. The electronic states correspond- 
ing to partly flat bands are the strongly localized states 
near the zigzag edges, called "edge states". Analytical 
properties will be discussed in the next section. 

In the case of q > N, the Landau levels are not per- 
fectly formed and the band structure at <f> = is almost 
unchanged, because the edges interrupt the cyclotron mo- 
tion of the electron (Fig. M (b)). For q < N, where the 
ribbon width is sufficiently wide compared with the cy- 
clotron radius, the Landau levels are nearly developed 
(Fig. @(c)). 

As an example for q <C N and higher commcnsurates, 
we show the energy band structure of (j) = | in Fig. |^ 
(d). The 2 x 4-Landau subbands are formed. Between 
the Landau subbands, we can see the additional disper- 
sion. The-states of these dispersions are also localized at 
the edgelI3li-3ll3, which originate from the cyclotron mo- 
tion of electrons and not for topological reason. It should 
be noted that the partly flat bands are formed at E = 
again. 




71-71 

k k 
FIG. 8. The energy band structures of armchair ribbon 

with N = 50 for (a) = 0, (b) = ;-L , (c) = j-L 

and (d) <f> ~ \ f° r = \- The dashed line in (b) is the energy 

structure at <f> = for comparison. 
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Next, we show the energy band structures of the arm- 
chair ribbon with N = 50 for 4> — 0, 555, j^q, \ in Fig. || 
(a)-(d). For <fi — 0, the profile of band structure is almost 
identical to the one of the graphite sheet ( Fig. || (a)). 
Here, we can not see partly flat bands at E = 0. 

In the case of q > N, the Landau levels are not per- 
fectly formed and the band structure at = does not 
change. For q < N, where the ribbon width is sufficiently 
wide compared with the cyclotron radius, the Landau 
levels are again almost formed as shown in Fig. |^. For 
the case of q <C N and higher commensurate fluxes, we 
show the energy bands structure of <f> = i in Fig. || (d). 
We also find that despite the additional dispersion be- 
tween Landau subbands, no partly flat bands, as may 
occur in zigzag ribbons, are present at E = even in a 
magnetic field. 



III. EDGE STATE 



The states corresponding to the partly flat bands are 
analytically derived by Fujita and coworker-s-for the semi- 
infinite graphite sheet with a zigzag edgeoB. It can be 
understood as localized states near the zigzag edge. It is 
also possible to find the edge states by solving the Harper 
equation fl2.Sp. 



At the beginning, let us rewrite the Eq.(2.£) to the 
transfer matrix form, 



where b m is b m — 1 for m = 1, N and b n 
us take open boundary conditions, i.e 



(3.1) 



*0 



for others. Let 



(3.2) 



and we impose the condition e = 0. One can obtain, 



=D 



k ' 



(3.3) 



where D k is — 2cos(|). Then the convergence condition 
\Dk\ < 1 is required, because otherwise the wave function 
would diverge in the semi- infinite graphite sheet. This 
convergence condition defines the region |tt < k < it, 
where the partly flat bands exist. The charge density is 
shown in Fig. ^ at (a) k — it, (b) 71", (c) it and (d) |tt. At 
k = it, the charge is perfectly localized at the zigzag edge. 
When the wave number deviates from k = tt, the electron 
gradually penetrates towards the inner sites. Finally, the 
electron states completely extend at k = §7r. 




FIG. 9. The charge density of the edge state at (a) k = n, 
(b) k = |7r, (c) k = g7r and (d) k = |-7r, where the radius of 
the circle means the magnitude of the charge density. 

Similarly, the case of finite magnetic field, using 
Eqs.(3.1) and (3.2) under the condition of e = 0, we can 



derive the wave function on the 11-th zigzag line as 
y n =K? =1 D k (i), 



where 

D k (i) -- 
For the rational flux 



_2cos(|-^±i 
1 



TT(f>) (i > 2) 

(i=l) 



(3.4) 
(3.5) 



= -, there is a relation of 
. q 

Dk{i + q) = Dk (i), so that in the case of semi-infinite 
graphite Eq.(3.4) might be rewritten as 



(A fc ) 



n-l 



where 



A fe = nUD k (i) 



(3.6) 



(3.7) 



Thus, the edge states are modified in the presence of a 
magnetic field. The condition of the convergence of wave 
function becomes < 1, which then defines the region 
of the flat band. Note that Afc has q internal degrees 
of freedom. In other words, there are q solutions which 
give the same value of Afc. For example, in the case of 
4> = j, there are 4 wave numbers which give Afc = 0, 
i.e, k = ±§7r, ±j7r, the charge density corresponding to 
each wave number is depicted in Fig. [l^. The charge 
density does not penetrate to the inner sites farther than 
up to the (4+l)-th zigzag chain, and there are 4 kinds of 
localized modes corresponding to the 4 internal degrees 
of freedom of Afc. Thus, the edge states in a magnetic 
field behave as the zero-field edge states with q internal 
degrees of freedom. 



G 



(a) (b) 




(a) (b) 




FIG. 10. The charge density of the edge state at (a) 
k — — §71", (b) k — —iir, (c) k — \-k and (d) k = |7r, where 
the radius of the circle denotes the magnitude of the charge 
density. 

Next, we discuss the DOS of the edge state in the ab- 
sence of a magnetic field, which will be used in the cal- 
culation of the Pauli susceptibility of zigzag ribbons in 
the later section. As we have seen in this section, the 
edge state penetrate to inner sites when the wave num- 
ber changes from ir to |-7r. If we consider the graphite 
ribbons with width N, two edge states which come from 
both side of edge will overlap with each other and develop 
the bonding and anti-bonding configurations. Since the 
magnitude of the overlap becomes larger when the wave 
number approaches |-7r, the band gap between the bond- 
ing and anti-bonding state formed by the two edge states 
gets larger toward k — §7r. Therefore the partly flat 
bands have a slight dispersion which depends on the rib- 
bon width N. In order to calculate the DOS, we first 
have to derive the precise energy dispersion for the edge 
states. The energy dispersion is calculated by the over- 
lapping of two edge states. From Eq.( |3.3| ), the amplitude 
of the edge state which penetrates from the first zigzag 
line is given by 

<S> n = D^- 1 = V Ai (3.8) 

which is located only on the A-sublattice. On the other 
hand, the amplitude of the edge state which penetrate 
from N — th zigzag line, is given by 

^N-n = D%~ 1 = *B, (3-9) 

which is located only on the B-sublattice. By using the 
tight binding Hamiltonian, the overlapping of two edge 
states is easily calculated, 



(y A \H\y B ) = ND£~ l (-2* - tD k ) = T k , (3.10) 

where D% = — 2cos(-|). Therefore, the energy spectrum 
of the edge states is given by the following eigenvalue 
problem. 

Uo*)(SM§) "■"> 

By diagonalization of this Hamiltonian matrix, we find 
the energy spectrum 

E k = -2tND% - 1 f-2t + 2t cos f J . (3.12) 

From this equation, around k — tt, the spectrum is given 
by E ~ k N . 

Therefore, the DOS related to the edge states has the 
form, 

8k 1 

P(e) = g~ e ~ ^, (3.13) 

where a — — 1. Note that this DOS has a power- 
law dependence, which is different from the ordinary van 
Hove singularity of p ~ observed in one-dimensional 
system. It is also found that the renormalized DOS is in- 
versely proportional to the ribbon width, which has been 
already confirmed by numerical calculation. □ 

IV. ORBITAL DIAMAGNETISM OF GRAPHITE 
RIBBONS 

It is well known that graphite shows a large anisotropic 
diamagnetic susceptibility, while aromatic molecules 
show only weak diamagnetism. This fact tells us that 
the orbital diamagnetic susceptibility is sensitive to the 
size of graphite fragments. Therefore, we would like to 
clarify the size and edge shape effect on the orbital sus- 
ceptibility, in order to understand the magnetic proper- 
ties of nanographite ribbons. Similar calculations of the 
orbital diamagnetic susceptibility. qn .carbon nanotubcs 
have been done by several authorsOO. They found that 
there are universal scaling rules in the orbital suscepti- 
bility as function of the Fermi energy, the temperature 
and the size of the nanotubes. Since it is also expected 
that graphite ribbons have such scaling rules, we study 
the scaling properties of the orbital susceptibility and the 
dependence on the edge shape. 

In this section, the orbital diamagnetic susceptibility 
Xorb of graphite ribbons is calculated in terms of the 2nd 
derivative of the free energy F(H, T) with respect to the 
magnetic field. The free energy is given by 

F(H, T) = fiN — -|- J &k^2ln(l + e-^- {H) -^ 

(4.1) 
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where (3 — and \x is the chemical potential and 

£k,n(H) (n is the band index) is the energy spectrum 
of the graphite ribbons in the magnetic field, as calcu- 
lated in the previous section. Then the magnetic moment 
M(H) and magnetic susceptibility x(H) P er s it e f° r finite 
temperature and arbitrary magnetic field H become 



^-km 



dk 



de 



•6 -4 



£-~> e /3(efc,n-M) 



1 dH ' 
(4.2) 



+ zigzag 
<•> armchair 







20 



W/a 



40 



and 



X(H) 



1 dM 
N~ P 'dH 



dk 



4 cosh 2 ^ 



de k , n 
dH 



1 



3/3(6fc,„-H) + l dH' 2 



(4.3) 



FIG. 11. The ribbon width dependence of the orbital dia- 
magnetic susceptibility \ort of graphite ribbons at T=0. 

The difference of the slope in Xorb and M or b is related 
to the more microscopic and configurational aspect of 
the ribbons. An important point to take into account is 
the ring current susceptibility in the equilibrium state, 
because the magnetic moments and ring currents are re- 
lated in the following way 



M = - I dVr ■ j. 



(4.5) 



The zero-field magnetic susceptibility Xo(T) for finite 
temperature is given by 



where r is the position and j is the current operator. 
J dV means volume integral. The current j is described 
on each bond by 



X(T) 



'iVeTT Go) e^--">+l 



dH 2 



Jij = i 



(4.6) 



(4.4) 



where S is the area of a hexagonal ring. N e is the electron 
number in the system. 

The width dependence of the orbital susceptibility Xorb 
at T = is shown in Fig. Analo gous to the graphite 
sheets, the aromatic molecules or the carbon nanotubes, 
the graphite ribbons exhibit diamagnetism. The magni- 
tude of Xorb{T rj 0) grows linearly with increasing W 
in accordance with the fact that Xorb of the graphite 
sheet diverges in the zero-temperature limit. A remark- 
able point is the different slope in the VF-dependence of 
Xorb for armchair and zigzag ribbons. Actually the dif- 
ference between the susceptibilities of the two types of 
ribbons increases for larger W. At a first glance, this 
result may seem unphysical, because one might attribute 
the difference to an edge effect which should diminish for 
wider ribbons. The origin of this discrepancy, however, is 
based on topological properties as we will show shortly. 
We would like also to mention the aspect that Xorb for 
the armchair ribbons shows some oscillations as a func- 
tion of W. This is due to the fact that armchair ribbons 
are metallic or insulating depending on W . On the other 
hand, no oscillations occur for zigzag ribbons, as they are 
metallic for all W . 



The ring currents also contribute to the linear response 
for weak magnetic fields. It is straightforward to cal- 
culated them and the corresponding susceptibility. The 
pattern of the ring currents for (a) a zigzag ribbon of 
N=10, and armchair ribbons of (b) N=18, (c) N=19 and 
(d) N=20 are shown in Fig. [l2](a)-(d), respectively. The 
magnitude (in the units of Jq(= 4t)) and directions of 
the ring current near the graphite edge is further given 
in detail in Fig 13. We can easily find that the current 
flow is symmetric with respect to x = and the total 
current in y-direction vanishes, because we consider an 
equilibrium state. The pattern of the current flow is strik- 
ingly different for the zigzag and the armchair ribbons. 
In zigzag ribbons, due to the lattice symmetry, the cur- 
rents along the vertical bonds are exactly zero. The cur- 
rents flow only along the horizontal bonds, whose direc- 
tions are antisymmetric with respect to x — 0. For arm- 
chair ribbons, the current distribution is quite different, 
because currents flow also on vertically oriented bonds, 
and exhibits a clear Kekule-type of pattern. Note that 
the Kekule pattern is more pronounced when N^3M- 
1(M=1,2,3,- • •), which corresponds to the semiconduct- 
ing armchair ribbon, whereas it is less distinct for the 
metallic case, N = 3M — 1(M=1,2,3,- • •), where almost 
no currents are found close to the ribbon center. The 
reason can be attributed to the interference effect be- 
tween the current flows assosiated with the two edges. 
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As shown in Fig. |l3](b), the currents are stronger along 
the cis-polyacetylene which is separated by one dimer 
lines. For N 3M— 1 and N = 3M— 1, the ring current 
patterns are depicted schematically in Fig. |lj (a) and 
(b), where thick bold lines denote dominant right-going 
currents while thick shaded lines are for the left-going 
currents. It is easy to find that for N ^ 3M — 1 both 
types of lines tend to avoid each other and form a Kekule 
pattern around the center of the ribbons. However, when 
N = 3M— 1, the lines lie perfectly on top of each other so 
that the left- and right-going currents cancel each other 
around the center of ribbon. Thus the effect of the lattice 
topology near graphite edge drastically changes the ring 
current flow in the whole sample. 



(a) 



(b) 



0.85 




2.5 



1.8 
0.2 
1.1 

Vj/~ 1 - 3 



M.3/ 

\0,2 



0.62 



FIG. 13. Current flow and magnitude near the edge for (a) 
zigzag ribbon (N=50) and armchair ribbon (N=50). 
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FIG. 12. The texture of the ring currents for (a) zigzag 
ribbon (N=10) and armchair ribbons of (b) N=18, (c) N=19 
and (d) N=20. In zigzag ribbon, because of the symmetry of 
the lattice, the ring currents along the vertical bonds are zero. 
In armchair ribbons of N=18 and 19, the Kekule pattern is 
clear. 



(a) 






FIG. 14. Schematic picture of ring current flow generated 
by the interference. 



(c) 




x/W u x/W " " D x/W 

FIG. 15. The position dependence of magnitude of the 
ring currents for (a) zigzag ribbons, (b) armchair ribbons 
(TV / 3M - 1) and (c) armchair ribbons (N = 3M - 1). 



In Fig. 15 



the magnitude of the ring current suscep- 
tibility J is plotted as a function of x/W with a fixed 
position y for (a) zigzag ribbons, (b) armchair ribbons 
with (TV 7^ 3M — 1) and (c) armchair ribbons with 
(N = 3M — 1). Interestingly, each graphite ribbon has a 
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scaling behavior as a function of x/W and the magnitude 
of the ring currents has a power-law decay. These facts 
further emphasize that the edge shape effect is significant 
in nanographitcs. 

Next we show the Fermi energy dependence of Xorb- 
Actually in real graphite materials, a small change in the 
carrier density from the half-filling is possible and can 
even be controlled by substrate properties. The calcu- 
lated Fermi energy dependence is shown in Fig. |l^, where 
it is found that Xorb/W is a universal function of (J.W. 
We normalize Xorb by dividing it by W, since it propor- 
tional to W ( Fig. |ll|). Furthermore we multiply Ep by 
W, because the direct gap at k = is proportional to the 
1/W at k — for armchair ribbons and at k = 2ir/3 for 
zigzag ribbons as is demonstrated in the Appendix. 





il 




, t * " A 


\ armchair(N=3m-i) J 


\ '■/ armchair (N=3m) .■' f 








B-B-H-H-B--H-H-H-B-H-H-B-H-H-H-* . 


• • • v " "•""»"•««"»»»>•«* 
x armchair(N=3m+i) 



-2 2 

EpW (units of t-a) 

FIG. 16. The Fermi energy dependence of the orbital mag- 
netic moments x°rb of graphite ribbons at T=0. 



As a final point in this section, we show the temper- 
ature dependence of Xorb in Fig. [I?], which is important 
from the viewpoint of experiments on nanographites. In 
all cases the magnitude of Xorb decreases with increasing 
temperature. It is also found that the temperature de- 
pendence of Xorb/W scales as a function of j3W, because 
the energy gap is proportional to the 1/W. Our calcula- 
tion also demonstrates that the edge effect becomes more 
significant with lower temperature. 




-armchair(N=3m-i) 
-armchair(N=3m) 
• armchair(N=3m+i) 



1 2 3 

W/(3(eV-a) 

FIG. 17. The temperature dependence of Xorb, where Xorb 
is scaled by 1/W and /3 is scaled by W . 



V. PAULI PARAMAGNETISM OF GRAPHITE 
RIBBONS 

In previous section, we have seen that the orbital dia- 
magnetic susceptibility depends on the edge shape in 
nanographite ribbons, especially, the topology of the lat- 
tice strongly affects the flow of diamagnetic ring currents. 
Here we discuss another important component of the 
magnetic susceptibility, Pauli paramagnetic susceptibil- 
ity XPi because zigzag ribbons have a sharp peak of DOS 
at the Fermi level. The width of the peak of DOS at the 
Fermi level has the order of meV, which is comparable 
to the temperature scale of room temperature. There- 
fore, it is expected that the Pauli susceptibility of zigzag 
ribbons might be sensitive to temperature, although the 
Pauli susceptibility of other usual metals is temperature 
independent. On the other hand, since the DOS of arm- 
chair ribbons at e = is zero or very tiny, we can neglect 
the effect of the Pauli paramagnetism in armchair rib- 
bons. 

The magnetic moment by Zeeman effect is 



AI = n B (n T - ni) , 



(5.1) 



where /is is Bohr magneton and rif (rij) means the elec- 
tron density with up-spin (down-spin). The electron den- 
sity at arbitrary temperature for each spin is given by 



1 



dfc 



1 



77 JlstBZ 



1 + cP( e n.k-°~fJ-B H) ' 



(5.2) 



where cr(=| , J.) means spin index. Therefore, the Pauli 
susceptibility xp per site is given by 



,. dM 
Xp = hm — — = 
A t^o dH irN e ^ 



dk 



1 



cosh(/3e n)fe )' 



(5.3) 



where f3 = -^f- Room temperature (T ~ 300K) corre- 
sponds to P ~ 0.25. We numerically calculated the finite 
temperature Pauli susceptibility of graphite ribbons us- 
ing this equation up to room temperature. 



20r 
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N=10 ° 06 
, N=20 

N=30 -\ 
.N=40 



N=50 



10 20 

1/|3 d0 3 eV) 

FIG. 18. The temperature dependence of xp f° r N= 10, 
20, .. , 50 up to room temperature. In the inset, the exponet 
is plotted on the width N and show a good agreement with 
the line of i — 1. 
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It is possible to clarify the contribution of the edge 
states to xp- As we have seen in Sec.HI, the DOS due 
to the edge stat es is given by E q.(3.15 ). After the sub- 
stitution of Eq.( 3.13 ) into Eq. (|5.3|) , we replace the k- 
integration by the energy integration. Then we can ob- 
tain the xp contribution due to the edge states as, 



1 



XP 



N e N(3° 



dx- 



cosh x + 1 



N ' 



(5.4) 



where x is (3ek and a is i — 1. Interestingly, \p nas the 
Curie-like temperature dependence, although in normal 
metals xp is basically constant in the temperature. The 
exponent of xp depends on the ribbon width through 
a. When N becomes infinite, the exponent a approches 
— 1 and xp show the Curie-law. However, in this limit, 
the con tribution of xp is diminished by a factor 1/N in 
Eq.@. 

Numerical results of the Pauli susceptibility xp °f 
zigzag ribbons up to room temperature are shown in 



Fig. 18 for various values of N. As expected, because 



of the edge states, xp shows Curie-like temperature de- 
pendence. In the inset of Fig.[l8], we plotted the N de- 
pendence of a, which was calculated by the least square 
method and has a good agreement with the line of — 1. 

The observed susceptibility x is essentially the sum of 
the orbital Xorb. and the Pauli susceptibility xp- The 
temperature dependence of the total susceptibility x is 
shown in FigjT^. The total susceptibility x shows the 
diamagnetic behavior in the high temperature regime 
and paramagnetic behavior in the low temperature. In 
the inset, the width dependence of the crossing temper- 
ature, i.e. x — 0, is plotted, which is well fitted by 

1 = 10 0.1526 x jj-LMfi 

Here we should remind that both aromatic molecules 
and bulk graphite show diamagnetic behavior, howerver, 
nanographite with zigzag edges have a remarkable para- 
magnetic behavior because of the edge state. If this para- 
magnetic behavior is experimentally detected, it will be 
an indirect evidence of the existence of the edge state. 



20 
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X N=40 
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FIG. 19. The temperature dependence of total susceptibil- 
ity x, which is Xorb. + Xp is shown for N = 10, 20, 50. 
In the inset, the width dependence of crossing temperature, 
where x — 0. 



VI. SUMMARY AND DISCUSSION 



In this manuscript, we discussed the electronic and 
magnetic properties of nanographite in magnetic field, by 
using the single orbital tight binding model with Peierls 
phase. Deriving the Harper equation for the graphite lat- 
tice, we studied the energy spectrum and dispersion in a 
magnetic field. At <j> = 0, it is found that zigzag ribbons 
have partly flat bands at the Fermi level, which is at- 
tributed to the edge localized states having non-bonding 
character. In terms of the Harper equation, the edge 
state can be analytically described even in presence of a 
magnetic field. We also studied orbital diamagntic prop- 
erties of the nanographite ribbons, where we found that 
the diamagnetic susceptibility Xorb. is very sensitive to 
the size and edge shapes of graphite ribbons. It is em- 
phasized that the flow of the orbital diamagnetic ring 
currents significantly depend on the lattice topology near 
the graphite edges. Especially, in the case of armchair 
ribbons, the pattern of the ring currents has drastically 
changed because of the interference effect of two edges. It 
is also found that the orbital diamagnetic susceptibility 
Xorb. is scaled as a function of the temperature, Fermi 
energy and ribbon width. Because the edge states in- 
duce a sharp peak of DOS at the Fermi level, the Pauli 
paramagnetic susceptibility should be an important com- 
ponent in nanographite with zigzag edges. Therefore, in 
the last section the Pauli susceptibility in zigzag ribbons 
has been studied, where we found that the zigzag ribbons 
with nanometer size show the Curie like temperature de- 
pendence of the Pauli susceptiblity in contrast to usual 
metals. From this significant contribution of the Pauli 
susceptibility, it is found that the observed susceptibility 
x(= Xp + Xorb.) of zigzag ribbons show diamagnetic be- 
havior at high temperature and paramagnetic behavior 
at low temperature. 

Here we introduce an interesting experimental re- 
sults, which might be connected with our theoreti- 
cal results. Some graphite-related materials consisting 
of nanographites, e.g., activated carbon fibers (ACF), 
amorphous carbons, carbon blacks, defective carbon nan- 
otubes etc., show actually anomalous behaviors in the 
magnetic susceptibility. While bulk graphite has a large 
diamagnetic and anisotropic susceptibility, a certain type 
of ACF with a huge specific surface area (SSA) up to 
3000to 2 /<7 (believed to consist of an assembly of minute 
graphite fragments with a dimension of 20 Ax 20 A) ex- 
hibits an paramagnetic response at room temperature 
and a strong Curie-like behavior in low temperature!!!! 
This kind of anomalous behavior of the susceptibility is 
also observed in many amorphous carbons and defective 
carpet-rolled carbon nanotubesEj. Although zigzag and 
armchair edges coexist in real carbon material, this be- 
havior of magnetic susceptiblity is consistent with our 
results. Although the sample production of graphite- 
related materials has still insufficient influence on size 
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and edge shapes, recently there are some experimental 
attempts to synthesize and nanographite systems and to 
control the size and edge shapes. One is "graphitization" 
of diamond powder with grain sizes 40 — 50A. Another 
method to produce nanogcaphites is epitaxial growth on 
substrates with step edgescJ. Therefore, we expect that 
in near future the magnetic properties of nanographites 
will be observable so that the influence of the edge states 
on magnetic properties can be tested in a controlled way. 
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APPENDIX A: ENERGY GAP 



In this Appendix, we analytically show that both the 
direct gap A a at k = of sufficiently wide armchair rib- 
bons and the direct gap A z at k = of sufficiently 
wide of zigzag ribbons are inversely proportional to the 
width of graphite ribbon W. This result supports that 
the Xorb/W is scaled as a function of the fj,W (fi is chem- 
ical potential). 

First, we examine the energy gap A a at k = of arm- 
chair ribbons. It is easy to find that at the k = the 
Hamiltonian can be rewritten as 

N 2 

H = -t^2(^2(a]^a j+ i^ + h.c.) + a)^a jfl + h.c.}, 

3=1 i«=l 

(Al) 

which is equivalent to the tight binding model for the 
2-leg ladder lattice having N rungaJ. The site in- 
dices (j,l) and (j,2) correspond to the jA(B) and jB(A) 
sites, respectively, when j is even(odd). The eigen- 
values are evaluated as = —2tcosmr/(N + 1) ± t 
(n=l,2,. . . ,N). It should be noted that the system is 
metallic only when Y=3m-1, because e + and e~ become 
zero for n — m and 2m, respectively. Therefore, A a 

are for Y=3to-1, 2 (2t cosi^^n) - t^j for Y=3m and 

2 ^cos(^±i.7r) - tj for N=3m+l, respectively. After 

the elimination of N in terms of W — (N— 3)^+-\/3 and 
the Taylor expansion under the condition of l/W <C 1, 
we can obtain the following results. 



W - 
n 

W 



s/3 

2 



N 
N 

N 



3M- 1 
3M 

3M + 1 



(A2) 



Thus the A a is inversely proportional to the ribbon 
width. 

Similarly, we can obtain the energy gap A z at k = ^ 
of zigzag ribbons. The Hamiltonian of zigzag ribbons at 



k = 4p is rewritten as 



H-- 



2N 

; E 

i=l 



(at a; 



h.c), 



(A3) 



which is equivalent to the tight binding model for the 
one-dimensional lattice having 2N sites. The site index 
i corresponds to iA, if i is an odd number, and to iB 
,if i is an even number. The eigenvalues are evaluated as 
e = -2icosmr/(2Y+l) (n=l,2,. . . ,N). Therefore, A z is 

7r) ] . After the elimination of N in terms of 



(2t cos( 



N+l 



2N+1 



w 



/3N 



1 and Taylor expansion under the condition 



of l/W <C 1, we can obtain the following results. 



A^ 



7T 

w 



(A4) 



Thus the A a is also inversely proportional to the ribbon 
width. 
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